PROCEEDINGS 

OF SCIENCE 



Strange quark mass and Lambda parameter by the 
ALPHA collaboration 



o 

(N 
o 

Q 

oo 



Marina Marinkovic 

Humboldt Universitat zu Berlin, Institut filr Physik, Newtonstr. 15, 12489 Berlin, Germany 
E-mail: jnarina . marinkovic@physik . hu-berlin . de 



Stefan Schaefer 

CERN, Physics Department, 1211 Geneva 23, Switzerland 
E-mail: 



Stefan . schaef er@cern . ch 



Rainer Sommer, Francesco Virotta 

NIC, DESY, Platanenallee 6, 15738 Zeuthen, Germany 



E-mail: Rainer . Sommer@desy . de , Francesco . Virotta@desy . de 



> 



We determine /k for lattice QCD in the two flavor approximation with non-perturbatively 
improved Wilson fermions. The result is used to set the scale for dimensionful quantities in 
CLS/ALPHA simulations. To control its dependence on the light quark mass, two different 
strategies for the chiral extrapolation are applied. Combining /k an d the bare strange quark 
mass with non-perturbative renormalization factors and step scaling functions computed in the 
Schrodinger Functional, we determine the RGI strange quark mass and the Lambda parameter in 
units of /k- 



CERN-PH-TH-2 011-314 
SFB/CPP-11-78 
DESY 11-241 
HU-EP-11/60 



XXIX International Symposium on Lattice Field Theory 
July 10-16 2011 

Squaw Valley, Lake Tahoe, California 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Strange quark mass and Lambda parameter by the ALPHA collaboration 



Marina Marinkovic 



1. Introduction 

The determination of the fundamental parameters of the standard model has a long tradition 
in lattice QCD. In particular the quark masses and the scale parameter A can be determined from 
first principles. This study is a part of a long-term programme of the ALPHA collaboration of 
computing these parameters, using the Schrodinger functional strategy to overcome the multi-scale 
problem and keep the full control over the systematic errors. 

The new ingredient presented here is the scale setting using a physical quantity, the Kaon 
decay constant /k- With this scale we achieve a 5% error, employing two different strategies for 
the chiral extrapolation which agree within errorbars. This enables us to give physical values for 
the RGI values of the strange quark mass and A-parameter, in the setup with two dynamical flavors 
of light quarks. 

The differences to the previously published values forMf GI |[I]] and A[|2|], for^Vf = 2, come from 
the improved scale setting. In the old computation, we used values available from the literature^], 
where the scale was set with r$. More recent determinations of r^/a find somewhat different 
results[Q, B]. Additional improvement comes from lattices with smaller pion masses and finer 
lattice spacing than previously available, giving a better handle on systematic effects. They were 
generated by the ALPHA Collaboration and the CLS 1 effort. 

2. Action and algorithms 

Our study is based on ensembles generated with the Wilson plaquette gauge action together 
with Nf = 2 mass-degenerate flavors of 0(a) improved Wilson fermions. The simulations are using 
either M. Liischer's implementation of the DD-HMC algorithm[0J or our implementation of the 
MP-HMC algorithm^]. 

The list of ensembles used in the analysis is shown in Table ||. Lattice spacings are ranging 
from 0.05fm to 0.08fm and their precise determination will be presented in the following section. 
The ensembles cover a wide range of pion masses going down to 270MeV, whereas all lattice 
volumes satisfy the requirement m n L > 4 to keep finite volume effects under control. 

3. Scale setting with /k 

To determine the scale and match to experimental values we have to extrapolate decay con- 
stants to the physical quark masses. For this we use two variants based on chiral perturbation 
theory(ChPT). The first one employs SU(3) chiral perturbation theory with a quenched strange 
quark. The aim here is to minimize the chiral corrections by keeping the sum M +M S of the light 
quark mass and the strange quark mass approximately fixed. Chiral corrections are expected to be 
well behaved, since in this setup all Goldstone bosons have a mass of at most the physical kaon 
mass (500MeV). The second approach uses heavy meson chiral perturbation theory (HMChPT), 
expanding only in the light quark mass (M = (M u +M a ) /2). Whereas the first strategy is most use- 
ful for Nf = 2, the second one is equally well applicable in Nf = 2 with a quenched strange quark 
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Table 1: Nf = 2 ensembles used in the analysis. MDU is the number of molecular dynamics units of the parts 
of the run chains used in the analysis. The configurations are saved after every 4 MDUs. As an illustration 
of autocorrelations, we give the integrated autocorrelation time of the pion mass expressed in MDU. T exp is 
estimated from /3 = 5.3 and quenched scaling^ and R act is the fraction of active links [S, ^|]. 

and in the Nf = 2 + 1 theory, the only difference being the low energy constants. The difference 
between the two strategies in approaching the physical point is illustrated in Figure [|. 

In the setup described in the next two sections we use two quarks with hopping parameters 
fCi = K2 = k and two additional quenched quarks with hopping parameters K3 = K4. 
From these we build pseudoscalars, pions with 
mass m n from two quarks with K\ and \K%. The 
kaons we build from (it, K3). The physical point 
is defined by m^ jP hy S = 134.8MeV and m K p h ys = 
494.2MeV, the values in QCD with the electro- 
magnetic interaction being switched off[Q]. The 
two strategies differ in how K3 is chosen as a func- 
tion of k. 

In the following computations we have in- 
cluded the effect of the autocorrelations in the er- 
ror analysis in a very conservative way. Namely, 
for the estimation of the error we take into account 
the tail of the autocorrelation function |gjp. Thus, 
we are convinced that we have statistical errors 
fully under control. The examples of autocorrela- 
tion functions for f K and /k are shown in Figure 
|3[ They are computed following the procedures 
detailed in 0]. 

3.1 SU(3) Chiral Perturbation Theory (Strategy 1) 

In this approach we define the strange quark hopping parameter K3 through the dimensionless 



strategy 1 ■ 
strategy 2 ■ 



Figure 1: A sketch of the two approaches for 
chiral extrapolation to the physical point. Strat- 
egy 1 imposes the condition on the sum of the 
strange and light quark mass: M s +M = const. + 
0(M 2 ), where M = (M u +M d ) /2, while the sec- 
ond strategy keeps the strange quark mass M s 
constant during the extrapolation. 
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for each value k and the gauge cou- 
pling (j8), such that the l.h.s. of the 
equation remains equal to the constant 



R 



n K,ph> 



phvs 



fi 



-. Rather than a fixed strange 



quark mass, this corresponds to M s + 
M = const, to lowest order in the expan- 
sion in the quark masses and this is ex- 
pected to give a flat chiral extrapolation 
for /k- The condition ([OJ) determines 
a value of K3 = h(k) as a function of the 
sea quark hopping parameter and it can 
be obtained by interpolation. After the 
dependence of K3 on k is determined, 
it remains to extrapolate the decay con- 
stant afa(k,h(k)) to the physical point, 




Figure 2: Chiral extrapolation of the kaon decay constant 
in lattice units. The values of /k have been multiplied by 
1 , but the error shown does not take the statistical error of 
the renormalization constant into account. The extrapolated 
results coming from both strategies agree at the physical 
point. They are obtained by applying the global fit for all 
three values of gauge couplings j3 , imposing a cut y n < 0. 1 . 

defined by the dimensionless ratio 



<(M(ft)) _ <, P hys 

/t(M(ft)) /i jPhy s' 



In the last step we use the prediction of this functional form coming from SU(3) ChPT[ 11]: 



af K (k,h{K)) = «/K,lat[l+lK(^,3'K) + («4-^)(^-^) + 0(y 2 )], 



(3.2) 



(3-3) 



£k(;Wk) = -^log(^) - ^l«g(^ - 1) + ^logOn) + iy^og{y„), (3.4) 
where af^\ at is the value of the decay constant in lattice units and the variables y are defined as 



"4 (ft) 

S^fl(k) 
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yK 



K,phys 
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yn 
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(3.5) 



' K,phys J K.phys 

The described chiral extrapolation to the physical point is shown in Figure |2| Finally, the lattice 
spacings for each value of the gauge coupling can be obtained with 



lat 



/K.phys 

and its values, together with the errors of this determination, are shown in Table [ 



(3-6) 



3.2 Heavy Meson Chiral Perturbation Theory (Strategy 2) 

In this approach, we work at the fixed strange quark mass and perform a chiral extrapolation 



in the light quark mass (M) using HMChPT [12]. Since we are working with Wilson fermions, one 
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Figure 3: Autocorrelation functions of pion decay constant(left) and kaon decay constant(right). In the error 
computation the tail of the autocorrelation function is taken into account. Standard truncations of the integral 
lead to an error smaller by a factor 2 or more. 



has to keep fixed the axial Ward identity (PCAC) mass of the strange quark. If we denote 



mi 



2 (PV pfi) ' 



(3.7) 



with P lj = \jfj(x)y$\ifj(x) and Aq T (x) = Yi( x )YaY5¥j{ x ) + a c hdoP'i{x) being pseudoscalar density 
and improved axial current, then the bare PCAC masses of the sea and the valence quark are defined 

by 

m\(k) = mi2 m 3 (K-,K- 3 ) = m 34 . (3.8) 

We first interpolate in K3 and determine functions s{k,fx) such that the strange quark mass is kept 
fixed to m^k, s(k,[i)) = jU. For fixed /I we then perform a HMChPT extrapolation to the chiral 
limit denned by ml phy Jfl phys ^ using 



3 

af K (k,s(k,n)) = P(lx)[\--[y K log(y Jl )-y 7l log(y Jt )} + a H (y Jt -y K ) + 0(M 2 
a 2 m 2 K (k,s(k,n)) = G(m)[1 + «h(^-^) + C(M 2 )]. 

In the end, the scale is obtained by interpolation in jit to the physical strange quark mass 

,2 



(3.9) 
(3.10) 



a 



K' . 2(Mk) 
at 



in 



K.phys 



f K - J P( AtK ) 2 /I, phys 



(3.11) 



The values of the lattice spacings from this strategy are shown in Table |2| Comparing to the results 
of the first strategy, we find a very good agreement as demonstrated in Figure ^[ 



a[fm] 



Asyst.fi! 



Strategy 1 


5.2 


5.3 


5.5 


0.0750 


0.0655 


0.04847 


0.0024 


0.0010 


0.00048 


0.0013 


0.0011 


0.00079 



Strategy 2 


5.2 


5.3 


5.5 


0.0745 


0.0649 


0.04808 


0.0025 


0.0010 


0.00047 


0.0014 


0.0012 


0.00090 



Table 2: Lattice spacings from the first strategy obtained by applying SU(3) ChPT(left) and from the second 
strategy, based on HMChPT(right). Estimation of the systematical errors is preliminary. 
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Figure 4: Chiral extrapolations of the product of lattice kaon decay constant /k and matching scale L max , 
defined with g 2 (L msx ) — 4.484 (left) and bare strange quark mass in units of fx (right). The matching scale 
in the extrapolation of ?77 s (right) is defined with g 2 {L max ) = 

4.61 



4. Determination of A and m s 

In the Nf = 2 theory, low and high energy physics have been connected non-perturbatively 
by the ALPHA Collaboration, using an intermediate (Schrodinger functional) renormalization 
scheme^ §]. Here QCD is formulated in a finite box of spatial size L and temporal extent T. 
The fields are subject to Dirichlet boundary conditions in time and periodic in space, where the 
former provide an infrared cutoff to the modes of quarks and gluons. This allows to perform simu- 
lations at zero quark mass and thus use the SF as a mass-independent renormalization scheme. We 
additionally specify that T = L and then the renormalization conditions are naturally imposed at 
the scale jJ. = I /L. 

To calibrate the overall energy scale, one fixes a large enough value of the coupling g 2 (L max ) to 
be in the low-energy region and relates the associated distance, L max , to a non-perturbative, infinite- 
volume observable, in our case fx- The extrapolation of the combination /k^max to the continuum 
limit is shown in Figure ||(left). It has been performed for both strategies of scale determination and 
the results agree within the errorbars. Combining the continuum result (/k£max)cont = 0.318(14)(6) 
and the value of (Aj^§-L max ) from [§] we get the updated value of the MS A-parameter in two flavor 
QCD: 

A £= it 7 , (AMs L max)/K = 316(26)(17)MeV, (4.1) 

\J K^max j 

where the matching has been performed at the low energy scale 1/Lmax> defined with g 2 (L max ) = 
4.484 and /k is the experimental value of the kaon decay constant. 

Furthermore, we compute the strange quark mass. We base it on the PCAC mass /Ik from the 



second strategy for chiral extrapolation (cf. Sect.p.2[). The MS strange quark mass is given by 



(small corrections proportional to the quark masses in lattice units are also accounted for) 

m% (2GeV) = — ^ ^s(2GeV) = 10L4(4 2)(2 5)MeV (4 2) 
*ms^ > m SF (L max ) Z P (L max )/ K M JK 

where the first factor is taken over from |[T]], and the new continuum extrapolation of the sec- 
ond factor is shown in the right panel of the Figure |J Here, the matching scale is defined with 
g 2 {L max ) = 4.61. The conversion factor to MS scheme = 2GeV) = 0.7431 is computed at 

4-loop perturbation theory; all other factors are non-perturbative. 
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5. Summary and outlook 

An important step missing in our previous work on non-perturbative renormalization of two 
flavor QCD has been performed. The scale is set from a physical quantity /k- Two strategies 
of chiral extrapolation are used and both give comparable results. Finally, we presented the non- 
perturbative computation of the A parameter and the strange quark mass of two flavor QCD. The 
only perturbative input in the whole calculation is the 4-loop conversion factor from the RGI strange 
quark mass to MS strange quark mass, which is required to make contact with wide spread conven- 
tions. 
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